Superfluid-Insulator transitions of bosons on Kagome lattice at non-integer fillings 
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We study the superfluid-insulator transitions of bosons on the Kagome lattice at incommensurate 
filling factors / = 1/2 and 2/3 using a duality analysis. We find that at / = 1/2 the bosons will 
always be in a superfluid phase and demonstrate that the T3 symmetry of the dual (dice) lattice, 
which results in dynamic localization of vortices due to the Aharonov-Bohm caging effect, is at the 
heart of this phenomenon. In contrast, for / = 2/3, we find that the bosons exhibit a quantum 
phase transition between superfluid and translational symmetry broken Mott insulating phases. 
We discuss the possible broken symmetries of the Mott phase and elaborate the theory of such a 
transition. Finally we map the boson system to a XXZ spin model in a magnetic field and discuss 
the properties of this spin model using the obtained results. 

PACS numbers: 



I. INTRODUCTION 

The superfluid (SF) to Mott insulators (MI) transitions 
of strongly correlated lattice bosons systems, described 
by extended Hubbard models, in two spatial dimensions 
have recently received a great deal of theoretical inter- 
est. One of the reasons for this renewed attention is 
the possibility of experimental realization of such mod- 
els using cold atoms trapped in optical lattices^. How- 
ever, such transitions are also of interest from a purely 
theoretical point of view, since they provide us with a 
test bed for exploring the recently developed theoretical 
paradigm of non-Landau-Ginzburg- Wilson (LGW) phase 
transitions^. In the particular context of lattice bosons 
in two spatial dimensions, a general framework for such 
non-LGW transitions has been developed and applied to 
the case of square lattice^. An application of these ideas 
for triangular lattice has also been carried oul£. 

The typical paradigm of non-LGW transitions that has 
been proposed in the context of lattice bosons in Refs. 
H andi is the following. For non-integer rational fill- 
ings / = p/q of bosons per unit cell of the underlying 
lattice (p and q are integers), the theory of phase tran- 
sition from the superfluid to the Mott insulator state is 
described in terms of the vortices which are non-local 
topological excitations of the superfluid phase, living on 
the dual lattice ££ These vortices are not the order pa- 
rameters of either superfluid or Mott insulating phases 
in the usual LGW sense. Thus the theory of the above 
mentioned phase transitions are not described in terms 
of the order parameters on either side of the transition 
which is in contrast with the usual LGW paradigm of 
phase transitions. Also, as explicitly demonstrated in 
Ref. 0, although these vortices are excitations of a fea- 
tureless superfluid phase, they exhibit a quantum order 
which depends on the filling fraction /. It is shown that 
the vortex fields describing the transition form multiplets 
transforming under projective symmetry group (projec- 



tive representations of the space group of the underlying 
lattice). It is found that this property of the vortices nat- 
urally and necessarily predicts broken translational sym- 
metry of the Mott phase, where the vortices condense. 
Since this translational symmetry breaking is dependent 
on the symmetry group of the underlying lattice, geome- 
try of the lattice naturally plays a key role in determining 
the competing ordered states of the Mott phase and in 
the theory of quantum phase transition between the Mott 
and the superfluid phases. 

In this work, we apply the theoretical framework de- 
veloped in Refs. to bosons on Kagome lattice described 
by the extended Bose-Hubbard Hamiltonian 
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at boson fillings / = 1/2 and 2/3. Here t is the bo- 
son hopping amplitude between nearest neighbor sites, 
U is the on-site interaction, V denotes the strength of 
the nearest neighbor interaction between the bosons and 
(i is the chemical potential. 

The main motivation of this study is two fold. First, 
since the geometry of the underlying lattice plays a signif- 
icant role in determining the nature of the Mott phase, we 
expect that such phases in Kagome lattice would be dis- 
tinct from their square 4 or triangular 5 counterparts stud- 
ied so far. In particular, the dual of the Kagome lattice 
is the dice lattice which is known to have T3 symmetry 8 . 
It is well known that the particles in a magnetic field 
on such lattice experience a destructive Aharonov-Bohm 
interference effect at special values of the external mag- 
netic flux leading to dynamic localization of the particles. 
This phenomenon is termed as Aharonov-Bohm caging in 
Ref. Ill In the problem at hand, the vortices reside on the 
dual (dice) lattice and the boson filling / acts as the ef- 
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fective magnetic flux for these vortices. Consequently, 
we find that at filling / = 1/2, the vortices become lo- 
calized within the Aharonov-Bohm cages& and can never 
condense. As a result, the bosons always have a feature- 
less superfluid ground state. Such localization of vortices 
and consequently the absence of a Mott phase for bosons 
is a direct consequence of the geometry of the Kagome 
(or it's dual dice) lattice and is distinctly different from 
expected and previously studied behaviors of bosons on 
square or triangular lattices^ whose dual lattices do not 
have T 3 symmetry. In contrast, for / = 2/3, we find that 
there is a translational symmetry broken Mott phase and 
discuss the possible competing ordered states in the Mott 
phase based on the vortex theory at a mean-field level. 
We also address the question of quantum phase transi- 
tion from such an ordered state to the superfluid and 
write down an effective vortex field theory for describing 
such a transition. 

The second motivation for undertaking such a study 
comes from the interest in physics of XXZ models with 
ferromagnetic J x and antiferromagnetic J z interaction in 
a longitudinal magnetic field Bi 

Hxxz = —Jx ^ (S?Sj + SfSj) + J z ^ SfSj 
(ij) (ij) 

-BiJ2 S ? (2) 

i 

where J x > and J z > are the strengths of transverse 
and longitudinal nearest neighbor interactions. Such 
spin models on Kagome lattice have been widely stud- 
ied numericallj*2*i2,. Further, the large J z limit of this 
model (in the presence of an additional transverse mag- 
netic field) has also been studied beforeii. A couple of 
qualitative points emerge from these studies. First, in the 
absence of external field Bi , such models on Kagome lat- 
tice do not exhibit S z ordering for any values of J z /J x . 
This absence of ordering is an unique property of the 
Kagome lattice and it has been conjectured that the 
ground state is quantum disordered^ Second, for B\ ^ 
and net magnetizations m = (S z ) = ft/6, ft/3, the model 
exhibit a quantum phase transition between a featureless 
state with (S x ) / to a translational symmetry broken 
ordered state with finite to. Using a simple Holstein- 
Primakoff transformation which maps the Hxxz (Eq. [2J 
to hardcore Bose- Hubbard model Hi (Eq. ^} , we show 
that it is possible to understand both of these features 
analytically at least at a qualitative level. The absence 
S z ordering for Bi — (which corresponds to average 
boson filling / = 1/2) turns to be a natural consequence 
of Aharonov-Bohm caging phenomenon discussed earlier. 
Further, the results from the analysis of the Boson model 
at / = 2/3 can also be carried over to study the possible 
S z orderings of the spin model at net magnetization ft/3 
which allows us to make contact with recent numerical 
studies in Rcfs. land Hi 

The organization of the paper is as follows. In the next 
section, we map between the spin model Hxxz to an bo- 
son model and carry out a duality analysis of this boson 



model. The dual Lagrangian so obtained is analyzed in 
Sec. UTTI for both filling factors / = 1/2 and 2/3. This is 
followed by a discussion of the results in Sec. IIVI Some 
details of the calculations are presented in Appendices lAl 
Eland© 



II. DUALITY ANALYSIS 

To analyze the spin model Hxxz, our main strategy 
is to map it to a Boson model using the well-known 
Holstcin-Primakoff transformation 

st = {sry = b i s;=(b\bi-^ (3) 

Such a transformation maps Hxxz to H\ given by 

H t = -t^^+h.c.) 

(ij) 

+ V J2(n i -f)(n j -f) (4) 

(ij) 

with the hardcore constraint m < 1 on each site. The 
parameters of the hardcore boson model Hi are related 
to those of Hxxz as 

Jx = 2t, J Z = V, B l = zv(f-^J (5) 

where z is the boson coordination number. In the Mott 
phases, the average boson density is locked to some num- 
ber and this will be / in Eqs0]and|S] In this work, we 
will consider the superfluid-insulator transitions with the 
average boson density fixed across the transition. 

In what follows, we would carry out a duality analy- 
sis of the hardcore boson model. Such a duality analysis 
of the boson model is most easily done by imposing the 
hardcore constraint in Hi (Eq. |2J| by a strong on-site 
Hubbard term potential to obtain -ffboson (Eq. [y. I n 
the limit of strong U the qualitative nature of the phases 
of the model is expected to be the same as that of the 
hardcore model. Hence for rest of the work, we shall con- 
sider the boson model H\ )OSOn . Also, since we are going 
to carry out a duality analysis of this model, we shall not 
bother about precise relations between parameters of Hi 
and i/boson, but merely represent to be a chemical po- 
tential which forces a fixed filling fraction / of bosons as 
shown in previous work 4 . 

The dual representatio n of Hi, oson can be obtained in 
the same way as in Refs. I ll r l The details of the duality 
transformation are briefly sketched in Appendix 1X1 The 
dual theory turns out to be a theory of U(l) vortices, re- 
siding on the sites of the dual (dice) lattice shown in Fig. 
n coupled to a fictitious magnetic field which depends on 
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FIG. 1: Kagome lattice (thin black line) and the unit cell of 
its dual dice (thick blue lines) lattice. The three inequiv- 
alent sites A, B and C of the dice lattice are shown in 
the figure. The sides of the rhombus ABAC have length a 
and angles 7r/3. The basis vectors u = (3a/2, \/3a/2) and 
v = (3o/2, — v/3a/2) that span the dice lattice are shown as 
arrowed lines. Note that the A sides reside in the center of 
the hexagons of the Kagome lattice while the B and C sites 
reside in the center of the triangles. 



the boson filling /. The dual action is given by 

Sd = ^2 X! ( £ A^A A„^4 fcA ~ /<Vr) 2 

P 



b . ft 



(6) 



where tpb are the vortex field living on the site b of the 
dual (dice) lattice, Ab^ is the U(l) dual gauge field so 
that t rv \ts. v Ab\ — rii where rii is the physical boson 
density at site i, denotes sum over elementary rhom- 
bus of the dice lattice, denotes lattice derivative along 
/it = x, y, t, and / is the average boson density. The mag- 
netic field seen by the vortices is therefore the physical 
boson density. Note that the vortex action Sd is not self- 
dual to the boson action obtained from i?boson- There- 
fore we can not, in general, obtain a mapping between 
the parameters of the two actions, except for identify- 
ing € T1J \/S. v Ab\ as the physical boson density^. There- 
fore, in the remainder of the paper, we shall classify the 
phases of this action based on symmetry consideration 
and within the saddle point approximation as done in 
Ref . H 

The transition from a superfluid ({if>b) — 0) to a Mott 
insulating phase in Sd can be obtained by tuning the pa- 
rameter r. For r > 0, we are in the superfluid phase. 
Note that the saddle point of the gauge fields Ab^ in ac- 
tion corresponds to e TI ^\A^Ab\ = f, so that the magnetic 
field seen by the vortices is pinned to the average boson 
filling /. Now as we approach the phase transition point 
r = 0, the fluctuations about this saddle point ({ipb) = 0, 



Aby = fx) increase and ultimately destabilize the su- 
perfluid phase in favor of a Mott phase with (ipb) 0. 
Clearly, in the above scenario, the most important fluc- 
tuations of the vortex field tpb are the ones which has the 
lowest energy. This prompts us to detect the minima of 
the vortex spectrum by analyzing the kinetic energy term 
of the vortices 



Hi 



kinetii 



-y 



ipb + h.c 



(7) 



where the sum over a is carried out over the six-fold 
(for A sites in Fig. ^) or three-fold (for B and C sites) 
coordinated sites surrounding the site b. 

The analysis of -f/kinetic amounts to solving the Hofs- 
tadter problem on the dice lattice which has been carried 
out in Ref. 0. We shall briefly summarize the relevant 
part of that analysis. It is found that the secular equa- 
tion obtained from ifkinetic for boson filling factor / is 
given byJi 
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where A, B and C denote incquivalent sites of the 
dice lattice, 7 = 27r/ = 2tt<P/<Po for the boson fill- 
ing /, and (p is the dual flux through an elementary 
rhombus of the lattice. Here ip = (ipA,'<pB,'4 > c) is the 
vortex wavefunction which has been written as ipb = 
ip(x,y) = ip(x, k) exp(iKn) with k = \fZk y a/2 and in- 
teger n — 2y/(v / 3a), a is the lattice spacing of the dice 
lattice, and e is the energy in units of vortex fugacity y v . 
For obtaining the solutions corresponding to £ / 0, we 
eliminate for ipB and ipc from Eq. [21 to arrive at an one- 
dimensional equation for iPa(% = 3ma/2, k) = ^(m, k) 



2 ros —(m -1/2) + 
3t 

-2 ros —(to +1/2) + 



2 cos (37m + 2k) ipA (to, k) 
iJjA(m - 1, k) 

ip A (m+l,K) (9) 



where we have used the fact that with our choice of origin 
the A sites have x = 3ma/2 with integer m. For ratio- 
nal boson filling / = p/q, Eq. [5] closes after translation 
by q periods. However, for q = 3g' [g' £ integer], this 
periodicity is reduced to g— . 

A key feature of the wavefunction %j> that we would 
be using in the subsequent analysis is its transforma- 
tion property under all distinct space group operations 



of the dice lattice. Therefore we collect all such transfor- 
mations here for any general filling /. Such operations 
for the dice lattice, as seen from Fig. Q involve trans- 
lations T u and T v along vectors u = (3a/2,^/3a/2) and 
v = (3a/2, — \/3a/2), rotations 7rn/3 for any integer n, 
and reflection I x about the x axis. It is easy to check 
that these are the basic symmetry operations that must 
leave the Hamiltonian invariant and all other symmetry 
operations can be generated by their appropriate combi- 
nations. Following methods outlined in Ref. 0, one finds 
the following transformation properties for the wavefunc- 
tions: 

Ta=u, v ■■ 1>(x, y) -> ip(x -a x ,y- a y )uo a ^ 2 ^ (10) 
I x :^(x,y)^r(x,-y) (11) 
R n/3 : ^ A {x,y) - ip A ([x + V3y]/2, [y - y/3x]/2) 

y, U j[{y 2 - x ' 2 )/ i +^ x y/ 2 ] 

i> B (x, y) -> Vc (jx + V3y}/2, [y - V3x}/2j 

x J( 9 ! -x 2 )/4+^ 9 /2] 

^c{x,y) -> ip B {[x + V3y}/2,[y-V3 x]/2) 

x J(y 2 -* 2 )/ 4 +^y/ 2 } (12) 

where (x, y) are the coordinates of the dice lattice and 
w = exp(27ri/), Here tp may denote either iPa(x = 
3ma/2,y = y/3na/2) or ips an d "0c at appropriate posi- 
tions within the unit cell. Note that all other above men- 
tioned symmetry operations can be generated by combi- 
nation of the operations listed in Eqs.^HEl F° r exam- 
ple I y — I x Rtt — IxI\/3 an d under this transformation 
ipA(x,y) -> ip* A (-x,y) and ip B ,c(x,y) -> ip* CB (-x,y). 
Some details of the algebra leading to Eqs. ^J^] is 
sketched in Appendix IB1 

In what follows, we are going to analyze the vortex 
theory for f — 1/2 and / = 2/3 which turns out to be 
the simplest possible physically interesting fractions to 
analyze. 

III. PHASES FOR / = 1/2 AND 2/3 
A. f=l/2 

The key difference between square and triangular lat- 
tices analyzed previously^and the Kagome lattice stud- 
ied in this work, comes out when we study the filling 
fraction / = 1/2. From Eq. we find that for / = 1/2, 
7 = 7r, so that the entire spectrum collapses into three 
infinitely degenerate bands e = 0, i-s/6. This situation 
is in sharp contrast to the square or the triangular lat- 
tices where the vortex spectrum for rational boson fillings 
/ = p/q has a fixed number of well defined minima. It is 
possible that vortex- vortex interaction may make the de- 
generacy finite, but it is unlikely to be lifted completely 
by interaction 8 . Therefore the vortex band has zero or at 
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FIG. 2: Aharonov-Bohm cage for a vortex on site of the dice 
lattice. The probability amplitude for the particle to be at 3 
vanishes due to destructive interference between amplitudes 
from the paths — > 1 — > 3 and — > 2 — > 3. Consequently the 
vortex can never propagate beyond the rim of the cage (black 
circles) and is dynamically localized. Analogous cages can be 
constructed for vortices on sites B and C— . 



least extremely small bandwidth with no well-developed 
minima at any wavevector which implies that the vortices 
shall never condense^. Consequently, we do not expect 
to have a Mott insulating state at / = 1/2, as seen in 
previous Monte Carlo studies Hi on Kagome lattice^. 
Comparing the parameters of H\ and Hxxz (Eq. 01 , we 
find that such an absence of Mott phase also implies that 
XXZ models on Kagome lattice will have no S z ordering 
for Bi = 0. This result was also obtained by numerical 
study of the related transverse field Ising model on the 
Kagome latticed. 

Physically, the collapse of the vortex spectrum into 
three degenerate bands can be tied to localization of the 
vortex within the so-called Aharonov-Bohm cages as ex- 
plicitly demonstrated in Refs. |£J An example of such a 
cage is shown in Fig. [21 A vortex whose initial wavepacket 
is localized at the central (white) A site can never prop- 
agate beyond the black sites which form the border of 
the cage. This can be understood in terms of destruc- 
tive Aharonov-Bohm interference: The vortex has two 
paths — > 1 — > 3 and — > 2 — > 3 to reach the rim site 
3 from the starting site 0. The amplitudes from these 
paths destructively interfere for / = 1/2 to cancel each 
other and hence the vortex can never propagate beyond 
the rim of the cage. Similar cages can be constructed 
for vortices with initial wavepacket at sites B and C&. 
This dynamic localization of the vortex wavepackets in 
real space, demonstrated and termed as Aharonov-Bohm 
caging in Ref. naturally leads to spreading of the vor- 
tex wavepacket in momentum space. Hence the vortex 
spectrum has no well-defined minima at any momenta 
leading to the absence of vortex condensation. 

The question naturally arises regarding the nature of 
the boson ground state at / = 1/2 (or equivalently 
ground state of XXZ model at Bi = and J z ^> J x ). 
The dual vortex theory studied here can only predict the 
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absence of vortex condensation due to dynamic localiza- 
tion of the vortices leading to persistence of superfluidity 
(or S x ordering for XXZ model) at all parameter values. 
Beyond this, numerical studies are required to discern the 
precise nature of the state. Previous studies of the trans- 
verse Ising model, which is related to the XXZ model, 
in Ref. [ll] have conjectured this state to be a spin liq- 
uid state. However, the Monte Carlo studies of the XXZ 
model— did not find any sign of spin liquid state. 



B. f=2/3 

In contrast to / = 1/2, the physics for / = 2/3 (or 
/ = 1/3) turns out to be quite different. In this case, 
we do not have any dynamic localization effect and the 
vortex spectrum has well defined minima. Note that for 
/ = 2/3, Eq. [5]closes after one period and the magnetic 
Brillouin zone becomes identical to the dice lattice Bril- 
louin zone without a magnetic fieldS.. As a result, Eq. [5J 
admits a simple analytical solution 

e = — (6 — 2 cos (2k) + 4 cos (k) cos(fc m )) 1 ^ 2 
ipA(x,y) — exp (iraik x aj2 + nn) (13) 

where we have k € (0,7r) and k x £ (0,47r/3), and k m — 
3k x a/2. The spectrum has two minima in the magnetic 
Brillouin zone at (k x a,n) — (0, 7r/3) and (2ir/3, 2tt/3). 
The eigenfunctions ip = (ipA,ipB,ipc) j corresponding to 
these two minima can be obtained from Eqs. 1131 and |S] 
They are 



ipi — exp(z7rn/3)(c, — c, 0) 
2nin 



H'2 



exp 



(c,0,-c) 



(14) 



where c is a normalization constant. Thus the low energy 
properties of the vortex system can then be characterized 
in terms of the field <I>: 



4> 



^i(x)(^i(x,t) +^ 2 (x)<^ 2 (x, t) 



(15) 



where Lp\(2) are fields representing low energy fluctua- 
tions about the minima "01(2)- As shown in Ref. 0, the 
transition from the superfluid to the Mott phase of the 
bosons can be understood by constructing a low energy 
effective action in terms of the ip fields. 

To obtain the necessary low energy Landau-Ginzburg 
theory for the vortices, we first consider transformation 
of the vortex fields ip under the symmetries of the dice 
lattice. Using Eas. 11011111211131 and US one gets 



T u : (fx -> (fii cxp(-i7r/3) tp 2 - 


-> if2 exp(i7r/3) 


T v : ipi — > ipi exp(i7r/3) ip 2 — > 


ip 2 exp(-i7r/3) 


4 


: <Pi(2) -> <fi(2) 


h 


■ <Pl(2) ^2(1) 


#27r/3 


: fl{2) -* ^1(2) 


R-k/3 


: <Pl(2) -* <P2(l) 



Some details of the algebra leading to Eq. ^] is given in 
App. E] It is also instructive to write the two fields ipi 
and if2 as two components of a spinor field \ — (Vi) ¥2)- 
The transformation properties of \ can then be written 

as 

T u -X^ xexp(-ir z tt/3) 
T v ■ X -> xexp(ir z 7r/3) 

Iy-X^ T xX* 
R2tt/3 '■ X - * X 
Rtt/3 -X^TxX ( 17 ) 

where T X y ig are the usual Pauli matrices. Note that 
i? 2 7r/3 plays the role of identity here where R^/3 mixes 
the two components of the spinor field x- 

The simplest Landau-Ginzburg theory for the vortex 
fields which respects all the symmetries is 



L {2) L (4) ^(6) 

J-iq, I J-/ ? , I -1-/., 



E [i(^ 

a=l,2 

4 4) = «fl^il 4 + l 



eA M )tp a \ +r\<p a \' 



+ w|^l| 2 |<^2| 2 



= w (^^ 2 ) 3 + h.c. 



(18) 
(19) 

(20) 
(21) 



The above lagrangian density can also be written in terms 
of the x fields. In particular the sixth order term be- 



Ea=± (X*TaXY 



This turns out to be the 



(16) 



comes 

lowest-order term which does not commute with r z and 
breaks the U(l) symmetry associated with the relative 
phase of the bosons. A simple power counting shows 
that L v 6 ^ is marginal at tree level. Unfortunately, the rel- 
evance/irrelevance of such a term beyond the tree level, 
which is of crucial importance for the issue of deconfinc- 
ment at the quantum critical point, is not easy to deter- 
mine in the present context for two reasons. First, the 
standard e and large N expansions does not yield reliable 
results for 2 + ID field theories, at least for N = 2 which 
is the case of present interest. Second, while it is known 
from numerical studies that such a cubic anisotropy term 
is relevant for the single vortex species case (for XY tran- 
sitions in 2 + 1 D), the situation remains far less from 
clear for the two vortex casein. We have not been able 
to determine the relevance/irrelevance of this term in the 
present work. 

If it so turns out that is irrelevant, the situation 
here will be identical to that of bosons on square lattice at 
/ = 1/2. The relative phase of the vortices would emerge 
as a gapless low energy mode at the critical point. The 
quantum critical point would be deconfined and shall be 
accompanied by boson fractionalizatiorA On the other 
hand, if is relevant, the relative phase degree of the 
bosons will always remain gapped and there will be no 
deconfinement at the quantum critical point. Here there 
are two possibilities depending on the sign of u and v. If 



FIG. 3: Mean-field vortex state for v > with y>i ^ 0. The 
arrows denote the phase of the vortex wavefunctions $ at the 
sites of the dual lattice. Note that all the sites of the Kagome 
lattice see the same vortex environment and are therefore 
equivalent. Hence we do not expect a density wave at the 
mean-field level. 




FIG. 4: Mean-field vortex state for v < with |^i||y>a| = 1 
and Relative = 7r/3. The notations are same as in Fig. [3] The 
magnitudes of the vortex wavefunctions are listed on the left, 
and A, B and C refers to corresponding inequivalent sites of 
the dual (dice) lattice where the vortices reside. Here there 
are six inequivalent sites for the Kagome lattice which are 
listed in Fig. [S] 



u, v < and w > 0, the transition may become weakly 
first order whereas for u > 0, it remains second order 
(but without any deconfinement) . The distinction be- 
tween these scenarios can be made by looking at the or- 
der of the transition and the critical exponents in case 
the transition turns out to be second order—. 

The vortices condense for r < signifying the onset of 
the Mott state and now we list the possible Mott phases. 
Our strategy in doing so would to be list the sites on the 
Kagome lattice which sees identical vortex environments 
(within mean-field) and then obtain the possible ground 
states based on this classification for a fixed boson filling 
/ = 2/3. We would like to stress that these results are 
obtained within mean-field theory without taking into 
account quantum fluctuations. 



FIG. 5: Mean-field vortex state for v < with | 1 1 V2 1 = 1 
and S rc i a tivc = 27r/3. This state also has six inequivalent sites. 
Notations are same as in Fig. 2] 

First let us consider the case v > 0. Here only one 
of the vortex fields condense and consequently the Mott 
state; breaks the symmetry between D and C sites of the 
dual (dice) lattice. A plot of the mean-field vortex wave- 
function for (ifi) ^ is shown in Fig. El The direction of 
the arrows denote the phase of the vortex wavefunction 
at the sites of the dual lattice while the amplitude of the 
wavefunctions are listed in the figures. We find that for 
this state all the sites of the Kagome lattice see identical 
vortex environment and consequently, within mean-field 
theory, we do not expect density-wave ordering for the 
bosons. However, the up and down triangles have differ- 
ent vortex environments and one would therefore expect 
increased effective kinetic energy of the bosons around ei- 
ther the up or down triangles. Consequently one expects 
an equal superposition state of bosons in which there is 
equal amplitude of the bosons in the sites around the 
up (down) triangles^. 

For v < 0, both the vortex fields have non-zero am- 
plitude i.e. \tp±\ = | ^2 1 ^ 0, and the ground state does 
not distinguish between the B and the C sites of the dice 
lattice. The relative phase between the vortex fields is 
fixed by . 

^relative = 2lTp/3 for W < 

= 7r(2p+ l)/3 for w > (22) 

for integer p. The plot of the vortex mean-field wave- 
functions are shown in Figs. [5] for w > and Fig. [6] for 
w < 0. The corresponding inequivalent sites, for w > 0, 
are charted in Fig. [7| We find that there are six in- 
equivalent sites (defined as those seeing a different vor- 
tex environment) in the Kagome lattice. We now aim to 
construct different density wave patterns based on two 
rules: a) the boson filling must be / = 2/3 and b) the 
equivalent sites shown in Fig.[S]must have the same filling 
throughout the lattice^. 

Such possible states are sketched in Figs. [7| and [H] The 
state denoted by <&i is shown in Fig. [7] and has a 3 by 
3 ordering pattern. Here all the red and blue sites are 
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empty (or spin down ) while the black and the green 
sites (both open and closed circles) are occupied (spin 
up sites). The numbers in the center of the hexagons 
denote the sum of magnetization (in units of fi/2) of the 
sites surrounding the hexagons. For the state <I>i shown 
in Fig. these take values and 6. 

More complicated states $2 and $3 with 9 by 9 or- 
dering pattern are also possible. These are shown in 
Fig. [H] For the state $2 the bosons are localized in red, 
blue, green (closed circle) and black (open circle) sites 
whereas the green (open circle) and black (closed circle) 
sites are vacant. The net magnetization of hexagons for 
$2 takes values 0, 4 and 2 as shown in Fig. |H| The state 
$3 can be similarly obtained from $2 by interchanging 
the occupations of the black and green sites while leav- 
ing the red and the blue sites filled. This has the ef- 
fect of 2 <-> 4 for the magnetization of the hexagons. 
Interestingly, any linear combinations of $2 and $3, 
<£(#) = cos(6>)$2 + sin(#)$3 for any arbitrary mixing an- 
gle 8, is also a valid density wave ordered state. The most 
interesting among these states turn out to be the $(7r/4) 
which has an 3 by 3 ordering pattern. Such a state corre- 
sponds a superposition of filled and empty boson sites on 
the green and black (both empty and filled) sites whereas 
the red and blue sites are filled. Here the sum of boson 
fillings of the sites surrounding the hexagons, takes val- 
ues and 3 as can be inferred from Fig. [5] Notice that 
the two states 3>i and $(7r/4) constructed here has the 
same long-range ordering pattern. 

To distinguish between the states $1 and <£>(7r/4), con- 
sider the operator A/), = X^ehcxagon n ^ where (A4) gives 
the sum of boson fillings of the sites surrounding the 
hexagon with center b. The values of (JV a ) for differ- 
ent states has been shown in Figs. and |H1 We note that 
the distinction between the states $1 and <&(7r/4) can be 
made by computing the values of (J2 a -^' a -^' a + u ) w hich 
measures the correlation between J\f operators at sites a 
and a + u where u is the basis vector for the dice lattice 
shown in Fig. Deep inside the Mott phase such cor- 
relations should vanish for $1 while it will remain finite 
for $(7r/4). Alternatively, one can also compute the dis- 
tribution of the hexagons with different values of (Af a ) to 
achieve the same go&\m. Such a distribution, computed 
in Ref. U seems to be consistent with the state $(7r/4) 
obtained here. 



FIG. 6: A chart of the inequivalent sites for the mean-field 
vortex state for v < with \ipi \ — \if2 \ — 1 and (9 rc i at i vo = tt/3 
denoted by red, blue, green (open and closed) and black (open 
and closed) circles. Inequivalent sites are defined as those 
which see a different vortex environment in the surrounding 
dual lattice. The relevant vortex environments are sketched 
in Figs. 0] and 




FIG. 7: One possible density wave states $1 at / = 2/3 
(m = 1/3 in units of h/2) when the black and green (both 
solid and empty dots) dots have all filled (or spin up) sites 
while the blue and red sites are empty (spin down). The 
centers of the hexagons are marked with 6 and respectively 
and shows the sum of magnetization (in units of ft/2)of the 
sites surrounding the hexagons. This state shows a 3 by 3 
density wave pattern as can be seen in the figure. A similar 
state can be obtained for / = 1/3 by replacing the bosons 
with holes and vice versa. 



Recently exact diagonalization study of XXZ model 
on Kagome lattice has been carried out in Ref. [1(J Their 
results concluded the existence of 3 by 3 patterned RVB 
state for m = ft/6. We note that the RVB state proposed 
in Ref. 10 shall have identical long range ordering pattern 
and Ea^a+n) to the state $(7r/4) obtained in the 
vortex mean-field theory. Of course, a true RVB state, if 
it exists, is beyond the reach of the mean-field treatment 
of the vortex theory. Obtaining such states from a dual 
vortex analysis is left as an open issue in the present 
work. 



IV. DISCUSSION 

In this work, we have applied the dual vortex theory 
developed in Refs. to analyze the superfluid and Mott 
insulating phases of extended Bose-Hubbard model (Eq. 
0) (or equivalently XXZ model (Eq. |2J) on a Kagome 
lattice at boson filling / = 1/2,2/3. The dual the- 
ory developed explains the persistence of superfluidity 
in the bosonic model at / = 1/2 of arbitrary small val- 
ues of t/U, t/V seen in the recent Monte-Carlo studies^, 
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and shows that dynamic localization of vortices due 
to destructive Aharonov-Bohm interference, dubbed as 
"Aharonov-Bohm caging" in Ref. 0, is at the heart of 
this phenomenon. This also offers an explanation of the 
absence of S z ordering in the XXZ model at zero longitu- 
dinal magnetic field for arbitrarily large J z /J x as noted 
in earlier studies^. In contrast for / = 2/3, we find that 
there is a direct transition from the superfluid to a trans- 
lational symmetry broken Mott phase. We have derived a 
Landau-Ginzburg action in terms of the dual vortex and 
gauge fields to describe this transition. We have shown 
that the order and the universality class of the transition 
depends on relevance/irrelevance of the sixth order term 
in the effective action. In particular, if this term turns 
out to be irrelevant, the critical point is deconfined and 
is accompanied with boson fractionalization. We have 
also sketched, within saddle point approximation of the 
vortex action, the possible ordered Mott states that ex- 
hibit a 3 by 3 ordering pattern. The results obtained 
here are in qualitative agreement with earlier numerical 
studies 9,11 on related models. 

We thank R. Melko and S. Wessel for helpful discus- 
sions and collaborations on a related project. We are also 
grateful to T. Senthil and A. Vishwanath for numerous 
insightful comments. This work was supported by the 
NSERC of Canada, Canada Research Chair Program, the 
Canadian Institute for Advanced Research, and Korea 
Research Foundation Grant No. KRF-2005-070-C00044. 

Noted Added: While this manuscript is being prepared, 
we became aware of a related work by L. Jiang and J. 
Ye±&. 




FIG. 8: Density wave state $2 at / = 2/3. Here the hexagons 
with black and green (both solid and open circles) dots have 
net boson filling 3 (or magnetization in spin language) . The 
blue and red sites are filled (spin up). The centers of the 
hexagons are marked with 0, 2 and 4 respectively and re- 
flects the sum of magnetization of the sites surrounding the 
hexagons. This state shows a 9 by 9 density wave pattern (full 
pattern not shown in the figure). Another state $3 with the 
same periodicity is easily obtained from $2 by simply revers- 
ing the direction of the spins (or changing from being filled 
to being empty and vice versa) on all the green and the black 
sites (both open and solid black and green circles) while leav- 
ing the blue and the red sites unchanged. This has the effect 
of 2 <-> 4 for the magnetization of the hexagons. Note that 
any linear superposition $(#) = cos(6)3>2 + sin(#)<I>3 is also a 
possible ground state with (m) = 1/3. The state $(7r/4) has 
a periodicity of 3 by 3. 



APPENDIX A: DUALITY TRANSFORMATION 

In this section, we briefly sketch the duality analysis 
which leads to the dual action (Eq. EJ) . We start from 
the Bose-Hubbard model H] )OSOU (Eq. 1). First, we follow 
Ref.|H to obtain an effective rotor model from H^ oson 
given by 

#rotor = —t ^2 C0S \^a4>i J + V ^ flihj 
{ia) {ij) 



E 



u. ,„ ,/ 

—PMH + ~^ n i [Tli - 1) 



(Al) 



where the rotor phase (0) and number (n) operators sat- 
isfy the canonical commutation relation <fii,hj = iSij, 

A a denotes lattice derivative such that A a 4>i = 4>i+ a — 4>i 
and i,a — x,y runs over sites of the kagome lattice. 

Next, following standard procedures'^, we write 
down the partition function corresponding to if r otor in 
terms of path integrals over states at large number of 
intermediate time slices separated by width At. These 
intermediate states use basis of hi and <pi at alternate 
times. The kinetic energy term of H m tor which acts on 



eigenstates of <f>i can be evaluated as 

exp [iAr cos (A a 4>i)} 
-> ex P [~ J m/(2<Ar) - iJiaAafa] 



(A2) 



where in the last line we have used the standard Villain 
approximation^ and Ji a are integer fields living on the 
links of the Kagome lattice. Integrating, over the fields 4>i 
and defining the integer field = (n*, Ji X , Ji y ), we get 
the link-current representation of the rotor model (Eq. 

Jin i,H 

-A^J iT J jV ] l[S(AM (A3) 

(ij) 1 

where / = [ijU + 1/2, A = VAt, and we have rescaled 
the time interval At so that e 2 = <Ar = 1/U6r. 
The constraint of vanishing divergence of the currents 
(A M Ji M = 0) in the partition function Z comes from 
integrating out the phase fields cj). We now note that 
this constraint equation can be solved by trading the in- 
teger current fields in favor of gauge fields A^^ which 
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lives on the links of the dual(dice) lattice and satisfy 
Jin = £^v\&vA b \. Note that e r „ A A„v4 bA corresponds 
to the physical boson density rij = J lT ■ In terms of these 
fields the partition function becomes 



^exp 



"2^2 E ( e ^x^uA b x 



(A4) 



Here we have dropped the term proportional to A and 
have assumed that it's main role is to renormalize the co- 
efficient e 2 — . Next, we promote the integer- valued fields 
Abp, to real valued fields by using Poisson summation 
formula and soften the resulting integer constraint by in- 
troducing a fugacity term y v cos (27rAt^). Further, we 
make the U(l) gauge structure of the theory explicit by 
introducing the rotor fields 9b on sites of the dual lattice, 



and mapping 2irA 



2-KAbn - A„0 b . This yields 



Z = / DA J 2?(9exp 

-Vv E C0S ( 2nA ^ - &^b) 



(A5) 



We note that exp(i9b) corresponds to the vortex creation 
operator for the original bosons^. Finally, we trade off 
exp {iOb) in favor of a "soft-spin" boson field ipb as in Ref. 
to obtain the final form of the effective action 

z = JvaJ Veexp(-S d ), 

E ( € l™X A " A b\ ~ fS^rf 



s d 



-y v J2(^b +t ,e 2 " A ^yJb + h.c. 

b 

-E^i'+h^i 4 ) 



(A6) 



APPENDIX B: TRANSFORMATION OF * 

Here we present the details of derivation of the trans- 
formation properties of ip under rotation by 7r/3. The 
rest of the transformation properties can be derived in a 
similar way. 

First note that any general term in the Hamiltonian 
^kinetic (Eq. 0) can De written as 



E 



C (aW,a«)^( 



ai 2) ,a|, 2) J«xi> 



a 



(1) a (2)) 



Here a, j3 are site index which can take values A,B or C, 
a^ 2 ) denote a near neighbor site of a^, and the phase 

(2) 

factor 9 = ^ J fa A y dy is obtained with the gauge A = 

A y = Hx and / = 4>/4>o where 4> — Ha 2 ^/i/2 is the dual 
flux through an elementary rhombus of the dice lattice 8 
and 4>q is the flux quanta. For example a typical such 
term can be 

J2r A (aiV^c (aW + l/2,a( 1 ) + ^/2 

xe 2,i/( a W+l/4) (B3) 

Now consider a rotation by an angle 7r/3. After the ro- 
tation, a typical term in the Hamiltonian becomes 

e ^(«: (i) ; «; (i) )^(^ (2) '% (2) ) 

xe 2 mf (a^^)(a^^)/^ (B4) 

where the rotated coordinates are given in terms of the 
old coordinates by 



at, 



a x /2 — V3a y /2 
a y /2 + \fla x j2 



(B5) 



Now from the structure of the dice lattice from Fig. [3] 
we know that such a rotation interchanges B and C sites 
while transforming A sites onto themselves. Then com- 
paring terms IB II and IB4I we see that we need 

(B6) 



^(a' xl a' y )e~^<^ = ^(a x ,a y ) 



where (3 takes values A,C and B for (3' — A, B, C. There- 
fore we have 



(a'«) - a (a ( 2 )) = 0(a« a< 2 >) 



-0(a' (1) ,a' (2) ) 



(B7) 



Fortunately one has a relatively straightforward solution 

to Eg. IB7I Rewriting a^ and a^ 2 ) in terms of a'^ and 
a ,(2) 

a x = a' x /2 + V3a' y /2 
a y = a' y /2-V3a x /2 (B8) 
and using Eg. IB2I we get, after some algebra 



a(a' x ,a') = -2irif 



(B9) 



Using Eos. IB9l and IB6I one obtains 



il)p(a x ,ay) 



_ piot(a x ,a y ) 



^(a\,a" y ) (BIO) 



(Bl) 



where 

S(a«,a (2) ) 



where a" x = (a x + ^/3a y )/2 and a" y = (a y — \/3a x )/2. 
This is Eq. I12I of the main text where we have used 
a x = x and a y = y for notational brevity. All other 
transformation properties can be obtained in a similar 
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APPENDIX C: TRANSFORMATION OF 
VORTEX FIELDS ip 

Here we consider the transformation of 951.2. First let 
us concentrate on the translation operator T u . Under ac- 
tion of T u , $ -> $' = $(x-u x , y-u v ) exp (4mfu x y/ . 
Note that for / = 2/3, u x — 3a/2 and y — \/3an/2, the 
exponential factor become unity. Thus one gets 

$' = exp [iiv(n— l)/3] (c, -c,0)<pi 

+ exp [i7r(rn — 1) + 2iit(n — l)/3] (c, 0, —0)952 

(CI) 

Comparing Eg. IC II with Eas. 1141 and 1151 one gets 



9?i ~~ * Vi exp(— i7r/3) 
9?2 — > -<^2 exp(-27ri/3) 



(C2) 



Similarly we get the transformation properties of 95^2 
under T v and I x . 

The transformation under rotation 7r/3 is slightly more 
complicated. First we need to note that for x = 'ima/2 
and y = \[2>naj2, one has 



J(v 2 -* 2 )/4W3^/2] = exp[-37ri(n-m) 2 /4 + ' 



(C3) 



Next note that with our choice of origin all the A sites 
are such that m + n is an even integer and in such cases 



exp [— 37ri(n — rn) 2 /4 + nin 1 = exp [m(m — 3n)/2] 

(C4) 



for all integers m and n. Hence under a 7r/3 rotation, we 
have using Eq. ^| 



= exp [— 7ri(m + 3ra)/6] (c, 0, — c)tpi 

exp [— 3mm/2 + mn/Q] (c, — c, 0)952 

x exp [— 7tz(to — 3n)/2] 
= exp[i7r(m + 2n/3)](c, 0, — c) ipi 

+ exp[i7rn/3](c, — c, 0) if2 (C5) 



Comparing Ea. lC5l with Eas.lT4landlT31 one gets (fU2) — * 
V2(i)- ^ similar analysis yields 951(2) — ► Vi(2) under the 
operation i?27r/3- 
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